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Abstract 

Coarse-grained Langevin-type effective field equations are derived for classical 
systems of particles. These equations include the effects of thermal fluctuation and 
dissipation which may arise from coupling to an external bath, as in the Brownian 
motion of a single particle, or which may arise from statistical fluctuations in small 
parts of an isolated many-particle system, as in sound waves. These equations 
may provide some guidance for the analysis of mesoscopic or microscopic molecular 
systems, or for systems of hundreds to thousands of subatomic particles produced 
in high energy nuclear collisions. 
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1 Introduction 



The theoretical description of Brownian motion of a classical particle is well known and 
widely discussed; it is textbook material [Q. It is used in numerous practical applications, 
such as evaluating reaction rates at finite temperature in systems where thermal fluctua- 
tions are of vital importance 0, ||, ||, |]. However, for continuous classical systems thermal 
fluctuations are not so frequently discussed, and the actual evaluation of fluctuating forces 
and dissipative coefficients for dense, interacting, classical, continuous systems is not so 
well known. 

We were originally motivated to this study by the physics of nuclear collisions at high 
energy. In these collisions many subatomic particles are produced, mostly pions which 
are the main carrier of the nuclear force. The collision may be viewed as a miniature Big 
Bang where soon after impact a large amount of the initial translational energy is put 
into particle creation and entropy production. This system can be roughly characterized 
by a temperature. As time goes on, the system expands and cools. Eventually collisions 
become so infrequent that thermal equilibrium is lost and the particles free-stream to 
infinity and are detected. The properties of a system of pions, numbering in the hundreds 
or thousands, at high temperature may have very interesting properties. For example, 
pion fields transform under an internal symmetry group rather analogous to spins in a 
magnet. During the early stage of the expansion, when the temperature is high, the field 
may collectively point in a direction different than it does in the surrounding vacuum. 
This is referred to as a disoriented chiral condensate (DCC) 0. One would like to have 
coarse-grained field equations to describe the fluctuation and dissipation of DCC domains 
and their inevitable coalescence and evolution into the sourrounding vacuum . A good 
description is lacking. 

Our goal here is to develop some understanding of the DCC problem by considering a 
collection of classical particles undergoing Brownian motion and generalizations of such. 
This might appear to be a simple problem but it seems not to have been discussed in the 
literature. In particular there are subtle issues relating to the nature of what one considers 
the heat bath. For a single particle it is relatively straightforward, for a collection of 
particles it is not. In one limit, the particles of interest may each be coupled to an 
external heat bath but they may also interact with each other via forces that are more 
slowly varying than the ones operative between the particles and the heat bath. In 
this case Langevin equations of motion can be obtained for time scales that are short 
compared to the interparticle interaction times but long compared to the interaction 
times with the heat bath. In another limit, there is no external heat bath; one averages 
over a macroscopically small but microscopically large number of neighboring particles 
and seeks a Langevin equation to describe the motion of these subsets of particles over 
times long compared to the force fluctuation times between subsets. In the real world 
there may be a continuum of interesting problems lying between these two extremes. 

In the following sections we will analyze the two limiting cases outlined above. In 
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both cases we will do coarse-graining to obtain an effective field equation of motion, in 
one and three dimensions. 

2 External Heat Bath 

2.1 Recollection of simple Brownian motion 

This case is discussed in many textbooks, such as and often repeated in the literature, 
such as [^. For a particle of mass m connected to a heat bath and moving under the 
influence of external potential fields, for example a small object suspended in a fluid or gas 
in a gravitational field, and moving in one dimension, the generalized Langevin equation 
is 




— oo 



Here v{t) is the velocity of the particle at time t. The force due to external fields is labeled 
by G. The force due to the heat bath has been separated into two parts: F = F + F', 
where F' represents the rapidly varying, random, part of the force whose average value is 
zero and F represents the slowly varying part whose average value is not necessarily zero. 
The separation of these two components depends on the coarse-graining time chosen. 
To be a useful coarse-graining, this time must be large compared to the characteristic 
correlation time Tcor of the force, but must be small enough on an observational time scale 
to record the desired coarse-grained trajectory of the particle. For example, if the particle 
under consideration is immersed in a gas of much lighter particles, Tcor is the average time 
between collisions with the gas particles which define the heat bath. The friction kernel 
is 

Kis) = P{F'{0)F'is)) (2) 

where f3 is the inverse temperature of the heat bath (we use units with ks = 1 throughout) 
and the averaging is carried out with respect to the heat bath. This is referred to as a 
generalized Langevin equation; the normal Langevin equation is obtained in the approx- 
imation that K{s) = 2a6{s). This may be considered as the limit of the somewhat more 
general case where 

K{s) = exp (-|s|/rcor) . (3) 

''"cor 

Simple Brownian motion in more than one dimension is obvious. 
2.2 Coupled Brownian particles 

Consider a set of Brownian particles referred to as which can move in one dimension 
and which are in thermal interaction with a heat bath referred to as H; see Fig. ^ The 
thermal interaction is mediated by a fluctuating force F^{t) which has a mean period 
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Heat Bath H 



Figure 1: A series of Brownian particles, B^_i, B^, B^^i, interacting with a ther- 
mal heat bath H with random thermal forces Fi^{t) (thin lines) and with their nearest 
neighbours via conservative forces (heavy lines). 



of Tcor- This time, Tcot, characterizes the relaxation time needed for H to reestablish its 
equilibrium configuration if it is perturbed by some sudden small change. We assume 
that the conservative force connecting Brownian particles is slowly varying compared to 
Tcor; we will study the change of position and velocity of the Brownian particles assuming 
that changes are small in a time interval r > Tcor- 

Suppose that at time t all are in thermal equilibrium with H. Consider a set of 
macroscopically similar systems forming an ensemble. Then the ensemble average of the 
fluctuating forces vanishes at this moment. 

{FMo = (4) 

The subscript indicates that the average is taken in thermal equilibrium. Due to a change 
in the position or velocity of some B^ the force may change and our system, both 
and H, may deviate from thermal equilibrium. We have to evaluate this deviation. 

Apart from the fluctuating random force -F)j(t) there is another slowly varying, not 
fluctuating, force G^{t). The forces connecting the Brownian particles are of this type. 
We will be interested in flnding the corresponding slowly varying part of the velocity 
of each S^. Integrating the equation of motion, mVfj, — G^{t) + F^{t), for a short but 
macroscopic period of time r we get, in a particular microscopic configuration, 

t+T 

m Mt+r) - vM -G^r+ J F^{t')dt' . (5) 

t 

Here we have taken into consideration that can be considered constant during the short 
period of time r. The microscopic configurations are not under our control; we know only 
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the ensemble average of similarly prepared systems. Taking the ensemble average of both 
sides of the above equation of motion we obtain 

t+T 

m (t;^(t+r) - v^{t)) = G^r + ^ {F,,{t'))dt' . (6) 

t 

In general {F^(t')) is not independent of the motion of B^, otherwise it would be always 
the same as the mean value (F^(t))o = in thermal equilibrium. We will evaluate the 
changes in {F^(t')) following the lines of [|^ (sect. 15. 7). 

Let us consider the change during the time interval from t to t+r'. The velocity of 
changes from to v^{t+T'). The motion of this particle affects its environment. If 
t' is small enough the mean force (F^) changes but still depends on its earlier value at t. 
If r' exceeds T^or, the heat bath will reestablish its thermal equilibrium and will be found 
with equal probability in any of its Q accessible states. Since the energy of the i?^ change 
the energy of the heat bath changes as well. The total energy of the heat bath changes 
by 

AE{t') = ... + AE,^,{t') + AE,{t') + AE^+i(r') + (7) 

where AE^{t') is the energy given to the heat reservoir H by the particle B^. Since we 
consider small time increments only the local environment of is relevant. The number 
of other Brownian particles influencing the heat bath is not infinite but extends to some 
distance comparable to the mean free path or the spatial correlation length. The number 
of states available to H changes from fl{Eo) to fl{Eo + ... + AE^^i + AE^ + AE^j_^i + ...). 
This change of energy of the heat bath will influence the Brownian particles connected to 
it and modify the populations of microstates r in the ensemble corresponding to B^. We 
will use this information to estimate the change of the fluctuating thermal force acting 
on B^. 

The equilibrium probability of the occurrence of a given microstate r for B^ is propor- 
tional to the corresponding number of states available to the heat bath. Here we assume 
that the total number of microstates of all Brownian particles -B^-i, -B^, -B^+i, ... together 
is negligibly small compared to the number of microstates of the heat reservoir. Now we 
can compare the probability distribution over the microstates W^i of B^ at time t and r' 
later. 

Wr^^jtW) ^ Wr^jt+T') _ n{Eo + ... + AE^.i + AE^ + AE^+i + ...) 

Wr,{t) Q{E,) 

= exp[/3(... + A^^_i + AE^ + AE^+i + ...)] (8) 

Here (3 = dlnQ/dE is the inverse temperature of the heat bath, assumed to be constant. 
This means that if more energy is made available to the heat bath the probability to 
populate a particular microstate r of B^ increases. This is true even if that particular B^ 
takes energy from the heat bath as long as its neighbours add more. 
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Now we can estimate how population probabilities change with time. 

Wr^it+r') = W^l^exp[P{... + AE^_^ + AE^ + AE^+, + ...)] 

« V^/;)[l + /3(... + AE^_i + AE^ + AE^+i + ...)] (9) 

We can also evaluate how the ensemble average of the fluctuation force changes during 
this time interval. 

r 

~ E [1 + /^(- + + ^^M + + -Ml. 

r 

= {[l + f3i... + AE^_^ + AE, + AE,+, + ...)]F^)o (10) 

Since the ensemble average of the fluctuating force vanishes for the thermal equilibrium 
distribution W^'^ the above expression for the fluctuation force reduces to 

{F^{t+T')) = f3{{... + AE^_^ + AE^ + AE^+^ + ...)F^)o . (11) 

In this estimate the change of the energy of the heat bath is still undefined. This energy 
is, however, simply the negative of the work done by the fluctuating force on the Brownian 
particle B^,. 

t+T' 

AE, = - J F,{t") v,{t") dt" (12) 
t 

Inserting this expression into that for the fluctuating force we obtain 

t+r' 

{F^{t+T')) = -P I dt" {F,{t+T')F,{t"))o v,{t") . (13) 
t 

We use the summation convention where a repeated index is summed over. In this case, 
the sum on v runs over the neighbors of B^; it is cut off by the range of the forces and by 
the flnite time interval. There is no need to do an ensemble averaging over the velocities 
because they are much more slowly varying than the fluctuating forces. Thus we have 
estimated the expectation value of the force for an ensemble, weakly deviating from a 
thermal equilibrium, via expectation values obtained in thermal equilibrium. This is, of 
course, just linear response theory. 

Now the equation of motion for the Brownian particle B^ can be cast in the form 

t+T t'-t 

mMt+r) - v^{t)) = G^T-13 j dt' j ds {F^{t')F,{t'-s))o v,{t'-s) . (14) 

t 
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The ensemble averaging in the expression for the interparticle force may be dropped 
because the coordinates change much more slowly than the velocities, which are again 
much more slowly varying than the fluctuating forces. 

It is straightforward to perform the above derivation in 3 dimensions. The forces 
and velocities become 3-dimensional vectors and with spatial indices labeled by 
Roman letters and where the bold Greek indices now label the position. The work done 
by particle will contain the scalar product f^-F^ of its velocity and the random force 
acting on it. Eqn. (0) takes the form 

t + T t'-t 

m{v'^{t+T) - v^it)) = G'^T-f3 J dt' J ds {F;{t')Fl{t'-s))o vi{t'-s) . (15) 

t 

The friction kernel is proportional to the correlation function of the random fluctuating 
force. 

= /3{F;it)Flit-s))o (16) 

Here it is assumed that this function is time translation invariant, as is normally the case 
in thermal equilibrium. This assumption can be relaxed. We will discuss the properies of 
time translationally invariant correlation functions a little later. 

Now it is useful to separate the fluctuating force into two components. The first com- 
ponent is just the average value of the fluctuating force in the slightly out of equilibrium 
ensemble. 

(Flit)) = - I - I ds K%{s) viit'-s) (17) 

t 

It leads to a damping of the velocity. Note that the right side of this expression involves 
averages defined with respect to the unperturbed thermal ensemble. We note that the 
integrand is appreciable only when s -C r because the coarse graining time r was chosen 
to be much bigger than the correlation time of the random forces. The velocity is slowly 
varying and so it may be evaluated with t replacing t' . (This is just the first term in a 
Taylor series expansion: v{t' — s) = v{t — s) + [f — t)v{t — s) + ■ ■ ■) . The upper limit of 
the s integral may then be sent to infinity yielding the approximation 



{Flit)) = - j ds K%{s) viit-s) 





The nonlocality reflects the time delay between the motion of the particles and the re- 
sponding force. 

The second component is the most rapidly fluctuating part and is defined by its zero 
average with respect to the actual out of equilibrium ensemble. 

F';{t) ^ F^it) - {Flit)) (19) 
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The equation of motion including both the drift and the fluctuating parts of the velocity 

oo 

G'^-Jds K^^{s) viit-s) + F';{t) . (20) 



IS 

dv^it) 
m — ^ — 
dt 



In this analysis the stochastic forces of the external heat bath act on every Brownian 
particle the same way. The correlation between two forces, i^^j^(s), under normal cir- 
cumstances is expected to decrease with increasing distance between Brownian particles 

and Sjy, and with the time difference s. The summation over the neighbours v of -B^ 
extends to infinity in all directions; however, the contribution of more and more distant 
neighbours is expected to be rapidly decreasing. 

2.3 Continuum limit 

Rather than attempting to solve the equations of motion for a macroscopic number of 
particles it is oftentimes useful to approximate them by a continuous medium. This is 
effectively coarse graining, and is accurate so long as the length scales of interest are large 
enough and the time scales of interest long enough. With this in mind, let us replace the 
discrete particle labels /x and v with continuous position variables x and y. Displacement 
of the particles in the gas will be denoted by ?7*(x, t). The sum over particle index is 



replaced by an integral over position. Divide both sides of eq. (20) by the average volume 
Vq per particle. The mass density is p = tji/vq. Forces per particle then become forces per 
unit volume and are denoted by a lower case letter. The Langevin equation then becomes 
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pf,\-K,t)=g\-K,t)- J d'y J ds 0{y,s) fj^ - y,t - s) + . (21) 



Here 

k''{y, s) = /3(f (x, t)f (x -y,t- s))o (22) 

is the correlation function for the force densities of the heat bath. 

If the medium is isotropic then the friction kernel must have the tensorial structure 

k'^iy, s) = h{y, s)yY + k^^iv, s) (5^^ - y'y^) (23) 

where y = |y|, and ki^ and k^ are longitudinal and transverse correlation functions. When 
we are interested in lengths greater than those characterizing the friction kernel we can 
evaluate the velocity at position x and take it past the y integration. (More generally, 
T] would be expanded in a Taylor series about y = 0.) Carrying out the averaging over 
angles gives 

pr7(x, t) = g(x, t)- j ds T7(x, t-s) j d^y {^^k^iv, s) + -k^iy, s)) + f'*(x, t) . (24) 
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Usually it happens that correlations fall off exponentially. Then the friction kernel may 
be parametrized by correlation times, correlation lengths, and strengths. 

kh{y,s) = f3{fl)exp{-s/Ti^-y/Xj^) 

kT{y,s) = /3(/|)exp(-s/rT-|//AT) (25) 

When our interest is on times greater than those characterizing the friction kernel and 
when the actual displacement velocity is slowly varying on those scales we can take the 
velocity outside the s integration to obtain 

p r7(x, t) = g(x, t) - pjri{yi, t) + fx, t) , (26) 

where the damping constant is 7 = |7l + |7t with 

7L = ^ Jds J £y exp (-s/tl - v/Xl) = 47rr(3)^^^^^ (27) 

and with a similar expression for 7t. This is a manifestation of the fluctuation dissipation 
theorem. 

The interesting feature of this Langevin field equation is that the damping term is 
linear in the displacement velocity. This is the form which is normally used in phenomeno- 
logical settings. However, this is certainly not the most general Langevin field equation 
for the displacement, as may be observed from the truncation of the Taylor expansion in 
y, and as we shall see in the next section. 



3 Internal Heat Bath 

In this section we remove the external heat bath and group all the particles of the system 
into small subsystems labeled A^. Focussing our attention on one subsystem we can 
think of the all the remaining ones as constituting a heat bath. There will be forces 
acting among these subsystems as sketched in Figure 2. These forces can be separated 
into a part which is rapidly varying on the time scale of interest and a part which is slowly 
varying. 

First consider a one dimensional system. During a particular time interval of duration 
r there will be a net momentum transfer p^-i^n from A^^i to yl^ and a net momentum 
transfer p^^n+i from A/^ to A/^^i due to the random forces. The net force experienced by 
Afj, due to these is 

The averaged equation of motion is just like eq. (14) except that now the force F orig- 
inates in the intersubsystem interactions, not with an external heat bath. We make the 
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Figure 2: A scries of small subsystems, A^^-i, Af,, A^+i, interacting with each other 
via random thermal forces (thin lines). The random forces are classified into nearest 
neighbour Fj^-^t), next nearest neighbour F^^)(i), ... , and so on. In addition each 
subsystem interacts with its nearest neighbours via conservative forces (heavy lines) too. 



assumption that there are no correlations between random momentum transfers involving 
different pairs of subsystems. That is, 



(p^,^+i(t)p,,,+i(t - s))o = Tr^K^'\s)5f,,, 



(29) 



where as before, K^-^\s) is expected to fall exponentially with s. Substitution into the 
equation of motion gives the following. 



miv 



t+T t'-t 

,,{t+T)-v^{t)) = G^T+ J dt' J ds (v^+i{t'-s) -2v^{t'-s) +v^_i{t'-s))K^'\s) . 

t 

(30) 

Here each subsystem has been assumed to have the same mass, m, for simplicity. 

Next we should allow for the possibility of random forces acting between next to 
nearest neighbors. The arguments exactly parallel those for nearest neighbors. The force 
is 



The correlation functions are 



(31) 



(32) 



This may be continued for next to next to nearest neighbors ad infinitum. The averaged 
equation of motion, taking all of them into account, is 

m{v^{t+T) - Vf,{t)) = Gu,T 

+ J dt' J dsJ2 {^t^+At'-s) - 2vi_,{t'-s) +v^_,{t'-s))K^''\s) . (33) 
t '^=^ 

The function K^'^^s) undoubtedly decreases rapidly with increasing a. 
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In taking the continuum limit we replace the subsystem labels // and a with position 
coordinates x and respectively. We divide the equation by the average length, of 
each subsystem. The sum over cr gets replaced by an integral over y. 

p {v{x, t+r) - v{x, t)) = g{x, t)r 

t+T t'-t OO 

+ j dt' J ds J dy k{y,s) (v{x+y,t'-s) - 2v{x,t'-s) + v{x-y,t'-s)) . (34) 

t 

Here k{y, s) is the continuation of K^^\s)/Iq. Generally, the velocity will be more slowly 
varying with position than the correlation function characterizing the random forces. If 
that is so then the difference v{x+y, t'—s)—2v{x, t'—s)+v{x—y, t'—s) is well approximated 
by y^dlv{x,t'-s). Then 

t+T t'-t OO 

p {v{x,t+T) — v{x,t)) = g{x,t)T + / dt' / ds dlv{x,t'—s) / dyy'^k{y,s). (35) 



To obtain a Langevin equation we proceed as before. Divide through by r, replace 
t' with t in the argument of the velocity on the right side, let the upper limit on the s 
integration go to infinity, and denote the displacement by the variable r]{x,t). 

oo oo 

pTi{x, t) = g{x, t) + j ds dlf]{x, t-s) J dy y"^ k{y, s) + f{x, t) . (36) 



The rapidly fluctuating force per unit length, /', is the deviation from the average value 
in the perturbed system. It is constructed by taking the continuum limit of the sum of 
random forces acting on the subsystem minus the average of those forces, similar to eq. 
(19). The average is just the second term on the right side of the above equation. 

1 oo 

ri:Flr\t)-{fix,t))^f'{x,t) (37) 

In general it is difficult to find a simple closed expression for /'. One way is to generate 
it from knowledge of the correlation functions like eqs. (29) and (32). 

When our interest is on times and lengths greater than those characterizing the friction 
kernel and when the actual displacement velocity is slowly varying on those scales we can 
take the velocity outside the integration to obtain 

pii{x, t) = g{x, t) + 7* dlf]{x, t) + fix, t) (38) 

where 

oo oo 

Y^JdsJ dyy^k{y,s). (39) 
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The essential difference between tfiis equation and the one for an external heat bath is 
the second space derivative acting on the velocity, which is absent in the latter case. The 
origin is translational invariance. The external heat bath imposes a particular frame of 
reference, whereas the absence of an external heat bath means that the system must be 
invariant under boosts of constant velocity. 

To generalize to more than one spatial dimension is straightforward. We assume that 
random momentum transfers between different pairs of subsystems are uncorrelated and 
that the system is rotationally invariant, for simplicity. 

¥^,^+Mu,u+p{t' - s))o = Tt''5^,J^,,K^\s) (40) 
This can be used to calculate the force-force correlation function. 

<T,p 

The sum over cr runs from (Xi, (T2, era = 1 to 00, and similarly for p, in order not to double 
count the number of pairs. Then 

- P{Fl{t')Fl{t' - s))ovi{t' - .) = E Hut's) - 2vl{t'-s) + vl_^{t'-s)) k\J\s) . 

(42) 

This leads to an equation of motion similar to eq. (33) with spatial indices i and j in the 
appropriate places and with the scalars n and a replaced with location vectors /u and cr. 
We take the continuum limit in the usual way. 

t+T t'-t 00 

+ j dt' j ds J dyidy2dy3 k'^{y,s) (v\x+y,t'-s) - 2v^{x,t'-s) + v\x-y,t'-s)^ . 
t 

(43) 

Because of rotational symmetry the kernel k has the same structure as in eq. (23). Expand 
the velocity to second order in y and integrate over all directions of y. 

p{v\^,t+T)-v\i^,t))=g\^,t)T 

t+T t'-t 

^ml I dsV%\^,t'-s) Jd^y y^ H{y, s) + Ak^iy, s)) 
t 

t+T t'-t 

+ - J dt' J ds 9,V-v(x, t'-s)J d'y y2 (A;L(y, s) - kT{y, s)) (44) 
t 
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Here y = |y|. Converting this to a Langevin equation in the now famihar way we get 

1 °° 

pT)(x,i) =g(x,i) + — y dsV(V-r/(x,i-s)) j d^y y" (ki^{y, s) - kriy, s)) 



oo 

+ ^Jds V^r7(x, t-s) J d'y y' {K{y, s) + ^k^{y, s)) + f'(x, t) . (45) 



The rapidly fluctuating force f may be constructed along the same lines as in one dimen- 
sion. 

In the very low frequency limit the velocities can be evaluated at time t and taken 
past the s integral. The Langevin equation then reduces to the Navier-Stokes equation 

= -VP + 77 + (C + V (V • v) (46) 

with the addition of the rapidly fluctuating force. By comparison we can determine the 
shear (77) and bulk {() viscosities. 

1 °° 

V = / / "^'^ (^^T + A;l) 


1 °° 

C = T^l'^'l^'yy' - 2A;t) (47) 


The requirement that the viscosities must be nonnegative places a restriction on the rela- 
tive magnitudes of the moments of the longitudinal and transverse correlation functions. 



4 Comparison of Field Equations and Conclusion 

An isotropic gas or liquid can be fully described by a scalar density 4> = — V-ry char- 
acterizing the local relative compression of the matter. The Langevin equation for the 
cases of external and of internal heat baths can be expressed in terms of alone without 
any reference to a vector displacement. But first we should inquire about the effective 
Lagrangian describing this field in the absence of dissipation and fluctuation. 
A real scalar field has the Lagrangian density 

£o = ^0^-^(V0)^-y(0) (48) 

where V{(f)) is a potential, usually a polynomial. This Lagrangian density does not contain 
the fiuctuations nor any dissipative term. For example, if V = then this describes 
undamped sound waves with speed c. The equation of motion is 

4) = c^V^(j)-dV{(j))/d(l). (49) 
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Here we may identify V-g with dV{4')/d(t). The influence of fluctuations can be included 
by adding to the Lagrangian density 

£fl„e = -^/(x,t)0. (50) 

The fluctuating field h' is simply identified with — V-f'/p in both the external and internal 
heat bath systems. The nature of this stochastic field is determined by the microscopic 
dynamics of the system expressed in terms of its original degrees of freedom. The frictional 
term cannot, of course, be written as an additional term in the Lagrangian. It must be 
added by hand to the equation of motion. 

Including the dissipative and the fluctuating terms the Langevin field equation for an 
external heat bath is 

0(x, t) = c2V^0(x, t) - dV{(t))/d(t) - 7(/)(x, t) + 6'(x, t) (51) 

where 7 was discussed in section 2.3. For the internal heat bath 

0(x, t) = c2v2(/)(x, t) - dV{(l))/d(f) + 7*V'(/)(x, t) + 6'(x, t) (52) 

where 

00 

^* = Y^pJ'^^J y' 5) + 2fcT(|/, s)) . (53) 

This is obtained under the assumption that the variation in is minor during the corre- 
lation times of /cl and /ex- 
it is important to appreciate the difference between the external and internal heat 
baths. With an external heat bath there is a preferred frame of reference. The random 
forces couple individual Brownian particles to the heat bath; it is assumed that there 
are no random forces between Brownian particles. With an internal heat bath there is 
no preferred frame of reference. The random forces couple different subsystems. This 
difference is the origin of the Laplacian in eq. (52) which is what really distinquishes the 
two field equations. In fact, if we desired a more fine-grained description we should expand 
the velocities in eqs. (15) and (43) to higher order in spatial and temporal variations. 
The resulting field equation will involve dissipation of the form 

E r,,(v^)^"^(|y 0(x,t). (54) 

Translational invariance demands that Fij — 0. All the coefficients can be determined in 
terms of moments of the friction kernel. Generally we would expect that only the first 

few terms in the sum are required for a good description of a coarse-grained system. 

Of course, in realistic situations it is always a matter of judgement whether there is 
a clear case of external or internal heat bath. It is possible to have different components 
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in the system, where part of this system can be considered as internal while the rest 
is external. Usually the external heat bath is chosen to be dominant, but this is not 
necessarily the case and mixed cases may come up in realistic studies. 

Although it is not our intention to apply these results to any particular problems in 
this paper it is instructive to consider one example. Let V{(j)) = and look for plane wave 
solutions in the absence of significant fluctuations. Denoting the frequency of the wave 
by uj and the wave vector by q we obtain the dispersion relation 



for the internal heat bath. For a real wave vector the frequency, in each case, has an 
imaginary part resulting in damping of the wave. Note that for the external heat bath 
the wave becomes overdamped for wave vectors less than 7/2c. This is certainly not 
representative of sound waves in the atmosphere! For the internal heat bath the damping 
goes to zero as resulting in the absence of dissipation in this limit so that uj —>■ cq. This 
dispersion relation is the same as obtained for damped sound waves in a gas or liquid by 
solution of viscous hydrodynamics 

In summary, we have derived (in some sense) coarse-grained Langevin-type effective 
field equations based on the classical dynamics of systems of many particles. We con- 
sidered two extreme limits: one where the stochastic forces arose from coupling to an 
external heat bath, the other where the stochastic forces arose from statistical fluctua- 
tions in small parts of the full system. The obvious next steps are to consider quantum 
mechanical effects, and to apply these results to real atomic, molecular, or nuclear sys- 
tems, equations 
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